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ABSTRACT 

In relativistic, Poynting dominated outflows, acceleration and collimation are 
intimately connected. An important point is that the Lorentz force is nearly 
compensated by the electric force therefore the acceleration zone spans a large 
range of scales. We derived the asymptotic equations describing relativistic, ax- 
isymmetric MHD flows far beyond the light cylinder. These equations do not 
contain either intrinsic small scales (like the light cylinder radius) or terms that 
nearly cancel each other (like the electric and magnetic forces) therefore they 
could be easily solved numerically. They also suit well for qualitative analysis 
of the flow and in many cases, they could even be solved analytically or semi- 
analytically. We show that there are generally two collimation regimes. In the 
flrst regime, the residual of the hoop stress and the electric force is counterbal- 
anced by the pressure of the poloidal magnetic fleld so that at any distance from 
the source, the structure of the flow is the same as the structure of an appropriate 
cylindrical equilibrium conflguration. In the second regime, the pressure of the 
poloidal magnetic fleld is negligible small so that the flow could be conceived as 
composed from coaxial magnetic loops. In the two collimation regimes, the flow 
is accelerated in different ways. We study in detail the structure of jets conflned 
by the external pressure with a power law proflle. In particular, we obtained 
simple scalings for the extent of the acceleration zone, for the terminal Lorentz 
factor and for the collimation angle. 

Subject headings: (magnetohydrodynamics:) MHD - relativity - galaxies:jets - 
gamma rays: bursts 



1. Introduction 

Highly collimated, relativistic jets are observed in active galactic nuclei (AGNs), mi- 
croquasars and gamma-ray bursts (GRBs). According to the most popular model, these 
outflows are powered hydromagnetically. By analogy with pulsars, it is assumed that the 
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magnetosphere o f a rapidly rotating accreti on disk (lLovelacdll976l : lBlandfordlll976l ) and the 
black hole itself (IBlandford &: Znajekl 119771 ) opens into a relativistic wind that transfers the 



energy away in the form of the Poynting flux. A long debated question is how and where the 
electro-magnetic energy is transferred to the plasma. The Poynting flux could be transferred 
to the kinetic energy of the flow by gra dual acceleration however, the acceleratio n strongly 
depen ds on the geometry of the flow (jChiueh et al.l Il99ll : iBegelman fc Lil Il994j : IVlahakis 
20041 ) so that acceleration and coUimation are intimately connected. 



General theorems affirm (iHeyvaerts &: Norman]ll989l . l2003l : IChiueh et al.lll99ll : iBogovalov 



19951 ) that at the infinity, the fiow should coUimate to the rotational axis, a good fraction 
of the electromagnetic energy being converted into the kinetic energy. However, it has been 
found that without an e xternal confinement, the characteris t ic coUimation /acc e leration scale 



IS ex p onentially large f Eichlei 19931: Begelman fc Li 1994 : Tominiatsu 1994 : Beskin et al. 
19981 : iBogovalovl \l9m Ichiueh et al1ll998l : IBogovalov fc TsinganosI ll999 L That is why in 



pulsar winds, the Poynting fiux is conv erted into the pla sma energy predominantly via dis- 
sipation processes (see, e.g., review by iKirk et al.l (120071 )). On the other hand, relativistic 
jets are observed in the sources where interaction of the outfiows with the external medium 
could not be neglected. In accreting systems, the relativistic outfiows from the black hole 
and the internal part of the accretion disc could be confined by the (generally magnetized) 
wind from the outer parts of the disk. A widely accepted model of long-duration GRBs 
assumes that a relativistic jet from the collapsing core pushes its way through the stellar 
envelope. In all these cases the external pressure could be responsible for coUimation of 
Poynting dominated outfiows. Moreover, the fiow is efficiently accelerated in the coUimated 
outfiows so that a significant fraction of the Poynting fiux could be eventually converted 
into the plasma kinetic energ y. Note that non-magnetized lets could also be efficiently fo 



cused by an ambient inedium (jEichleiill982l : iPeter fc Eichler 



1995: Levinson fc Eichler 2000 



Bromberg fc LevinsonI 120071 ). An advantage of magnetically driven outfiows is a relatively 
low mass load, which naturally leads to highly relativistic velocities. 

An explicit solution for the relativisti c magnetiz e d wiri d from the accretion disk was 
found in the force- free approximation by IBlandford Jl976h. In this solution, the mag- 
netic surfaces are nested paraboloids. iBeskin &: Nokhrinal (120061 ) generalized this solution 
to include the inertia forces and showed that the magnetic surfaces are only slightly mod- 
ified and that the fiow is accelerated until the equipartition level. A few se U-similar solu- 



tions to the relativ i stic magnetohydrodynami c ( MHD) equations wer e found (ILi et al.lll992 



ContopoulosI I1995I : Iviahakis fc Konigll boOSal B | Narayan et al.l 120071 ) , which resemble out- 
fiows from a disk. These solutions also demonstrated that coUimation and acceleration 
could oc cur at a reasonable, even though large, scale. Nunaerical simulations support these 
findings jKomissarov et all 120071 . l2008l : lTchekhovskov etaP 120081 ). 
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A crucial assumption in these models is a non-zero magnetic flux threading the disk 
and the black hole. The total flux should in fact be inflnite (going to infinity with the outer 
disk radius) becaus e it is the pressure of the poloidal field, not the hoop stress, that colli- 
mates the outfiow (ISpruit et al.l 119971 ) . Such a field could not be generated in the disk; it 



should be dragged inward by the accreting material (IBisriovatvi-Kogan fc Ruzmaikirull976 



Bisnovatyi-Kogan fc Lovelacd 120071 : iRothstein fc Lovelacd l2008l ) . Magnetized outfiows with 
the ze ro net magnetic fiux, the so called magnetic towers, were proposed by iLvnden-Belll 
( 1996) and then studied b oth analytically (Lovelace fc Romanovall2003]: Uzdenskv fc MacFadven 



20061 : iLvnden-Belj 



2006 



Sherwin &: Lynden-Bell 



and numerically (jhovelace et ah 2002; Kato et al 



2007 



2004 



Gourgouhatos fc Lvnden- Belli boosh 



Nakamura et al 



200e 



d. l2007h 



Since 

there is no large scale magnetic field in this model, the jet is collimated by the pressure of 
the ambient medium so that an extended outfiow surrounding the jet is anyway needed. 

It is well known that in relativistic MHD outfiows, the acceleration zone spans a large 
range of scales. This is because the electric force, which is negligibly small in the non- 
relativistic case, becomes comparable with the Lorentz force and when the fiow velocity 
approaches the speed of light, these two forces nearly cancel each other so that both accel- 
eration and coUimation proceed very slowly. Within the light cylindei0, the magnetosphere 
corotates with the central source so that the plasma, which slides along the rotating field 
lines, could acquire only moderate relativistic velocities. Beyond the light cylinder, the fiow 
is accelerated at least until the velocity exceeds the fast magnetosonic velocity. The fast 
magnetosonic point is already very far from the light cylinder but in this point, the plasma 



energy is still well below the Poynting fiux (e.g.. iCamenzindl (119861 )). The complete transfor- 
mation of the electro-magnetic to the kinetic energy could occur only at the scale much larger 
than even the distance to the fast magnetosonic surface. This means that a few different 
spatial scales are present in the problem, which poses a strong challenge to numerical simu- 
lations. On the other hand, multi-scale systems are suitable for asymptotic analysis. In the 
spirit of the method of matched asymptotic expansions, one can solve the equations in two 
overlapping domains, namely, in the near zone, f2r ~ 1, where the force-free approximation 
could be used, and in the far zone, where one can considerably simplify the equations in the 
limit (Qr)^^ ^ 1. Both solutions should be matched in the intermediate region where the 
fiow is still force-free but the condition fir ^ 1 is already fulfilled. 

In this paper, we study properties of relativistic jets at the distances much larger than 
the light cylinder radius. First we obtain the asymptotic equations describing the fiow in 
the limit Qr ^ 1. Far enough from the source, these equations are valid till the axis of the 



^In differentially rotating magnetospheres, the surface fir = 1 is not a cylinder but we retain the standard 
term, which has come from the pulsar theory. 
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flow so that these equations in fact describe the whole flow in the far zone. We apply the 
obtained equations to jets confined by an ambient medium. We show that there are two 
different regimes of the fiow coUimation and acceleration. In the first regime, the structure 
of the flow at any distance from the source is the same as in an appropriate cylindrical jet, 
i.e., the residual between the magnetic hoop stress and the electric force is compensated by 
the pressure of the poloidal field. We will refer to this regime as to equilibrium coUimation 
in the sense that the fiow remains in the cylindrical equilibrium. In the second regime, one 
can neglect the pressure of the poloidal field so that the dynamics of the fiow is the same as 
in the case of purely toroidal field; this regime will be called non-equilibrium. In different 
coUimation regimes, the acceleration regimes are also different. 

We show that while the flow is Poynting dominated, the structure of the jet is governed 
by a simple ordinary differential equation, which could be easily solved for any distribution 
of the external pressure. The general theory will be applied to jets with a constant angular 
velocity propagating in a medium with the pressure decreasing as a power law. We also 
study the structure of the moderately magnetized core of the jet; such a core is presented 
near the axis of even Poynting dominated flows because the Poynting flux vanishes at the 
axis. As the jet propagates, the flow is accelerated and the inner parts of the jet reach 
equipartition between the kinetic and electromagnetic energy so that the the moderately 
magnetized core expands within the jet. Depending on the proflle of the conflning pressure, 
the core could either occupy only internal part of the jet so that the main body of the flow 
remains Poynting dominated or expand till the periphery of the flow such that the whole jet 
ceases to be Poynting dominated. 

The paper is organized as follows. In the next section, we shortly outline derivation of 
the basic equations describing relativistic, axisymmetric MHD flows. In Sect. 3, we shortly 
discuss the boundary conditions and integrals of motions. In Sect. 4, we flnd asymptotic 
equations for the flow in the far zone. In Sect. 5, we use the derived equations to develop 
a technique for flnding the structure of coUimated, Poynting dominated jets. In Sect. 6, we 
apply this technique to jets with a constant angular velocity propagating in a medium with 
the pressure decreasing as a power law. The terminal Lorentz factor of the flow as well the 
terminal coUimation angle, are estimated in Sect. 7. In Sect. 8, we study the structure of 
the moderately magnetized core of the jet. The obtained results are summarized in Sect. 9. 



2. Basic equations 



For the sake of consistency and in order to introduce nota t ions, let us shortly review the 
basic theory of relativistic, magnetized winds (jOkamotdll978l : iLovelace et al.lll986l : iLi et al. 
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19921 ). Let the plasma be cold, which is a good approximation in the far zone where the flow 



is already expanded. Then the steady state equation of motion is written as 

P7(v ■ V)7v = [(V ■ E)E + (V X B) X B] ; (1) 

where p is the plasma proper density, 7 the Lorentz factor, v the plasma velocity; the speed 
of light is taken to be unity. Here the second pair of Maxwell's equation is already used. 
The equation of motion should be supplemented by the flrst pair of Maxwell's equations, 

V-B = 0; VxE = 0; (2) 

by the continuity equation, 

V ■ (p7v) = 0; (3) 

and by the condition of flux freezing, 

E + V X B = 0. (4) 

In axisymmetric conflgurations, the magnetic fleld is conveniently decomposed into the 
poloidal and toroidal components, B = Bp + B^cj), the poloidal fleld being expressed via the 
flux function ^ 

Bp = -V^ X 0. (5) 

r 

Here cylindrical (r, 4>, z) coordinates are used; hat denotes unite vectors. The condition of 
flux freezing implies that the flux surfaces are equipotentials, which yields 

E = -^](^)V^; (6) 

where n(^) is the angular velocity of the fleld line. This gives a useful relation 

E = rVLBp. (7) 

The plasma streams along the flux surfaces so that the flow velocity may also be decomposed 
into the poloidal and toroidal components, v = Wpl + t>00, where 1 is the unit vector along 
the magnetic surface, 

T=fixJ; ii = V^/|V^|. (8) 

The condition of flux freezing yields a relation between the components of the velocity and 
magnetic fleld: 

BpV^ - B^Vp = rfi(^)5p; (9) 

which implies that the plasma slides along the rotating magnetic fleld lines. The continuity 
equation could be integrated, with the aid of Eq. ([2]), into the conservation law 

Anpvp-f = vWBp; (10) 
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where the function rj describes the distribution of the mass flux at the inlet of the flow. 



The three remaining equations are obtained by projecting the equation of motion onto 
directions 1, and n. The first two may be manipulated into the integrals of motion 



7 - 



V 



(11) 



(12) 



representing conservation of the energy and of the angular momentum, correspondingly. 
Note that the widely used parameter a, defined as the ratio of the Poynting to the matter 
energy flux, is presented via the basic quantities as 

/i-7 



a 



1 



(13) 



The projection of the equation of motion onto the normal to the flux surface, n, yields the 
transfield force-balance equation (the generalized Grad-Shafranov equation) 



1 

n 



2 2 I 
P7 + 



4:71 



Bl 



n ■ V— + — p7 v^n ■ r 



in 



;n.V[r\Bl-E'' 



(14) 



where IZ is the local curvature radius of the poloidal field line (defined such that TZ is positive 
when the flux surface is concave so that the coUimation angle decreases). 



1 = -fi-(T- v)T= ii-[rx((v xT))] = -0-(v xT). 



(15) 



Eqs. ([7]), ([9]), (fTOj) . (fTTj) . (fT2|) and form a complete set of equation describing cold, 
axisymmetric MHD flows. This set could be reduced to a pair of equations for \E' and 7. 

Eliminating B^ from Eqs. ffTTl) and (|T2l) . one can express the azimuthal velocity via ^ 
and 7 as 



1 



1 



7 



(16) 



Assuming for simplicity that at the origin of the outflow, the rotation velocity is well below 
the speed of light, f^rm, f (/,,in 1, one reduces Eq. (fT6l) to the form 

1 



^11-^ 



7 



(17) 



where the index "in" is referred to the parameters of the injected plasma. Substituting this 
relation into Eq. ([9]) and eliminating B^ with the aid of Eq. (|TT|) . one gets the expression 
for the poloidal velocity 



r'Q'Bp 
r]{fi - 7) 



1 - 



1 



1 - 



Tin 

7 



'1^ 
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Now one can write the identity + f ^ + 7 ^ = 1 as the BernouUi equation 



1 - 



Tin 

7 



+ 



1 - 



Tin 

7 



+ 



1; 



(19) 



which connects the Lorentz factor of the flow with the geometry of the flux tube defined by 
the function 

The transfield equation f[T^ is converted into an equation for \1/ and 7 upon substituting 
E from Eq. from Eq. ([TT]) . p from Eq. flTUl) . t>p from Eq. f[T5]) and f,^ from Eq. f[T7|l . 

Therefore Eqs. (1141) and (fT9!) form a complete set of equations. 



3. Boundary conditions and integrals of motion 



At the inlet of the flow, one should specify the distribution of the poloidal flux or, which 
is the same, of the poloidal magnetic field Bp. We are interested in outflows subtending 
a finite magnetic flux, \E'o) therefore we have to prescribe a boundary condition at the last 
magnetic surface. If the flow is confined by the pressure of the external medium, the pressure 
balance condition should be satisfied at the boundary. In the proper plasma frame, the 
magnetic field is B' = {B^ - E'^f'^ = {Bl + (1 - ^^2)52)^2. The condition that the 
pressure of this field is compensated by the external pressure is written as 



[Bl + (1 - r]V2)fi2]^(,,,)=*„ = 87rp,.t(r, z)- 



(20) 



where Pext is the pressure of the external medium. 



In the cold flow, one has also to prescribe the functions 0(^'), ri{^) and 7in(^) at the inlet 
of the flow so that only two integrals of motion, and /i(^'), remain unknown. Assuming 
that the rotation velocity is non-relativistic at the origin of the flow, we have eliminated the 
dependence on I (see transition from Eq. flTB]) to Eq. f[T71) ). In the general case, the integral 
I may be expressed via 7m(^) ^(^), vi^)^ /^(^) and Bp at the inlet of the flow making 
use of ([9]), (fTT!) and f|T2|) . So one has to find only the energy integral, /i. This integral is 
determined by the condition of the smo oth passage of the flow through the singular surfaces , 
Alfven and modified f ast magnetosonic (ILi et al. I I1992I : iTsinganos et al.lll996l : lBogovalovlll997l : 
Vlahakis et al.l l200d : IVlahakis fc Konigll l2003al ). In the Pointing dominated outflows, the 
Alfven surface coincides with the light cylinder, fir = 1, whereas the fast magnetosonic 
surface goes into the far zone fir ^ 1. Transition through the Alfven surface could be 
studied in the force-free approximation, i.e. neglecting the plasma energy and inertia. 



The force-free limit of the transfleld equation is obtained by taking p = in Eq. (I14p . 
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Making use of Eq. ([7]) , one finds 



1)^ + ^fi ■ V [{n^r^ - 1) Bl] = ■ V irB.r - 

In the force-free limit, the energy equation ( fTTi) is reduced to the form 

rB^ = 2I{^y, 21{^) = r]{^)fi{^)/n{^y, 



Q'Bln ■ r. 



(21) 



(22) 



which means that the current flows along the magnetic surfaces. Now the force-fr ee balance 
equat ion can be recast in the form of a second order elliptical equation for (jOkamoto 
19741 ). which is sometimes called the pulsar equation. By inspecting Eq. fl2T|) . one can see 
that in the pulsar equation, the second derivatives are multiplied by (fi^r^ — 1) so that the 
equation is singular at the light surface. The condition of regularity at th i s surf a ce enables 
one to fix th e poloidal curren t J(\l^) f e.g.. iFendt ( 1997 ): IContopoulos et al.l (119991 ): lUzdensky 
kooA bm^ : llovelace et"aD (120061 ): Iximokhinl ( |20od )). Then the energy integral is found 
just adding the matter energy fiux as /i = 7in + 2VLI /rj. The first term here is small in 
the Poynting dominated outfiows however, one cannot neglect it close to the axis where the 
current I goes to zero {I = tt J jrdr = (7r/2)j(r = 0)r^ — * as r — > 0). 

Note that decreasing of the energy fiux towards the axis is the generic property of the 
Poynting dominated outfiows because the poloidal current, /(^E'), always goes to zero at 
\1' ^ 0. Such a "hollow cone" energy distribution accounts, in particular, fo r a specific 



morph ology of the inner Crab and other pulsar wind nebulae (e.g., review by iKirk et al. 



( 120071 )). In any case, the exact shape of /^(^E') depends on the geometry of the fiow close to 
the source. In this paper, we study the fiow in the far zone therefore we assume that this 
function is given together with other integrals of motion. 

We would like only to note that the function /i(\I') has a universal form close to the axis. 
The poloidal field remains finite at the axis so that Eq.([2]) yields 

1 



^ = -r^Bjr = 0,2 



$ ^ 0. 



(23) 



Beyond the light surface, the magnetic field becomes predominantly toroidal whereas the 
fiow becomes predominantly poloidal (see the next section), therefore Eq. yields B^ ^ 
E = rVLBp. Then the second term in the energy equation is written, close to the 
rflB^/r] = Vt^r'^Bp{r = 0,z)/ri = 2[fi(0)]^\E'/?7(0). So close to the axis, the energy integral 
has the universal form 

/x(^) = 7i,(0) f 1 + i 1 : ^ ^ 0: 



where 



1 + 
7in(0)r/(0) 

mm' ■ 



(24) 
(25) 
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Note that the flow is Poynting dominated only at ^ \E'. 



An important point is that in outflows with a constant angular velocity, one can assume 
for the estimates that the energy integral is described by the linear function not only close 
to the axis but a c ross the jet. Both an anal ytical solution for the parabol oidal flux surfaces 



(iBlandford 



2007 



200 



1976 



Beskin fc Nokhrinal |2006| ) and numerical simulations (IKomissarov et al. 



Tchekhovskoy et al.ll2008l ) show that this is a good approximation for such jets. 



One can also obtain a quite general estimate for the energy integral taking into account 
that beyond the light cylinder, Eq. (Q yields 5^ ~ —QrBp, which simply means that each 
revolution of the source adds to the wind one more magnetic loop. Then the second term in 
Eq. ( fTTI) may be estimated as {Qr)'^Bp/r]. Making use of the estimate \E' ~ {l/2)r'^Bp (the 
coefficient is exact when the poloidal field is homogeneous), one finds finally 



/i(^) ~ 7in + 



(26) 



This expression provides a rough estimate for the energy integral for arbitrary f2(\l') and 

77(^). 



4. The basic equations in the limit ^ 1 

4.1. Expansion in 1/r. 

We are interested in outfiows initially dominated by the Poynting fiux. In such outfiows, 
the Alfvenic surface, where ^ Bp, nearly coincides with the light surface Qr = 1. In the 
far zone, i7r ^ 1, the toroidal field decreases as oc 1/r , see Eq. (ITTi) . The poloidal field 
decreases as 1/r^ therefore in the far zone, the field is nearly toroidal. The fiow in the far 
zone becomes nearly radial because according to Eq. f|T7|) . oc 1/r. In spite of this, one 
generally have to retain the terms with Bp and in the equations. The physical reason is 
that the hoop stress is nearly compensated by the electric force so that one cannot generally 
neglect small pressure of the poloidal field. The formal reason is that the leading order terms 
in Eqs. ( IT9|) and (fl^ are the same, which makes the system nearly degenerate, so that one 
have to retain smaller order terms. 

In the transfield equation 0141) . the leading order terms are those in the right-hand side 
because the terms in the left-hand side are small either as Bp/ B^ or as r /TZ. In the Bernoulli 
equation f|T9l) . one gets the leading order terms just neglecting the terms with 1/r and I/7. 
This yields 

//(/. - 7) = n\''Bp- (27) 
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or, according to Eqs. ([7]) and (|TT1) . 



Bs + E = 0. 



(28) 



If one substituted this relation into the right-hand side of the tran sfield equation, one would 
kill the leading order terms. The correct procedure (1Vlahakisll2004j ) is to expand the Bernoulli 
equation (|T9|) to the first non- vanishing order in 1/r and I/7 and only then to eliminate the 
leading order terms from Eq. f|T4|) . Expanding Eq. f|T9|) yields 



vif^ - 7) 

Qr 



^2 



1 1 b;. 



(29) 



Substituting this relation into the right-hand side of Eq. ([Hj), one gets 



1 

n 



2 2 

P7 f p 



^2 



52 

p 



3 
47r 



2 2 
P7 ^</, 



n- r 



rli- V 



7' 

/ in 



(30) 



In this equation, there are no terms which nearly cancel each other. Therefore one can now 
retain only terms of the lowest order in 1/r and I/7, e.g., neglecting Bp with respect to E 
or substituting Vp by unity. Moreover, one can now use Eq. fl28l) . which is the zeroth order 
approximation to the Bernoulli equation, in order to further simplify this equation. For 
example, the expression in the first brackets in the left-hand side could be transformed, with 
the aid of Eqs. (171). ([TUD and ([TTD. as 



2 2 I 
PI Vp + 



E^ 



__ = _(4vrp7 Vp - QrBpB^) = 



(31) 



The expression in the second brackets in the left-hand side of Eq.fl30l) could also be simplified 
in the same way after substituting from Eq. ([TT]) : 



o2 

P I 2 2 



1 



-nrBpB^ + ArrpYvp 1 - — 



7ir 



Br, 



-VLrB^ + ?77 



^ _ Tin 

7 



T]Bp 



7 



27in+ — 

7 



(32) 



We can also use Eq. (127]) in the right-hand side of this equation. Eventually one finds 



PVBp ^ r]Bp 



n 



7 



H - 27in + ^ 1 K • r = —a ■ V 



2r2 



^272 



fi2r^ 



(33) 



This is the asymptotic transfield equation valid at fir ^ 1. It may be significantly simplified 
in specific cases. 
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4.2. Asymptotic transfield equation in different regimes. 

If the flow is initially Poynting dominated, one can neglect the terms with far enough 
fro m the axis, n, f ir 3> 7in- Then one comes to the equation obtained, in a different form, 
bv lviahakisi fl2004h : 

The solution to this equation describes the main body of the flow but it could not be 
continued to the axis because it could not satisfy the condition \l'(r = 0) = 0. Close to the 
axis, where the Poynting flux decreases according to Eq. 1^^, the terms with should be 
retained therefore one should generally solve the full asymptotic equation (133]) . 



Taking into account that the terms with play role only close to the axis, where the 
angular velocity and the injection Lorentz factor could be considered as constants, one can 
present Eq. fl55]) , with the aid of Eq. flTTI) , in a more convenient form 



1 

n 



n ■ r 



1 



2% 



+ 



1 
2^ 



1 + 



7-' 

/ m 



n ■ V 



(35) 



The solution to this equation could be continued to the axis in spite of the fact that the 
equation was formally derived at the assumption Vlr ^ 1. The reason is that far enough 
from the center, the light cylinder fir = 1 is well within the matter dominated zone, <C ^, 
where the flow is practically hydrodynamic. 

In the most interesting case of collimated flows, z ^ r, one can take n ■ r = r and 
n ■ V = 8/ dr. When looking for the shape of the magnetic surfaces, one can conveniently 
use the unknown function r{^,z) instead of \E'(r, z). Then, e.g.. 



Br, 



IV^I 



19* 

r dr 



dr \ 



(36) 



In the same approximation, the curvature radius may be presented as (note that TZ is defined 
to be positive for concave surfaces) 

1 d^r 

Now the transfield equation for the collimated flows in the far zone could be written as 



r]ii 



+ 



1 



27ir 



+ 



2r 



1 + 



7-' 

/in 



d rj^{fi — 7)' 



^2^2 J Q^, ^2^2 



(38) 



We believe that this equation suits well to numerical solution because it does not contain 
terms that nearly cancel each other. In many cases it could even be solved analytically. 
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For analytical solution, this equation could be conveniently considered in two overlapping 
domains, namely, in the main body of the jet, where the flow is significantly accelerated so 
that one can neglect the terms with 7in/7, and close to the axis, where the flux surfaces 
are nearly straight so that one can neglect the curvature term d'^r/dz^. Solutions in these 
domains are smoothly matched in the intermediate zone. 

Close to the axis, where the flux surfaces are nearly straight, 

(fr I 

d^^Qh^ ^^^^ 
one can neglect the term with the d erivative in z and write the t ransfield equation as an or 



dinary differential equation (see also lBeskin fc MalyshkinI (120001 ) : iBeskin fc Nokhrinal (12006 



20081)) 



Hi.f)-.,.^^(.-.)^^. (40) 



We will analyze it in sect. 8. In some cases, the condition (|39|) is fulfilled across the whole jet; 
then the full jet structure is described by the one-dimensional equation, the z dependence 
entering only via the boundary conditions. 

Note that neglecting the derivative in z in the the transfield equation, one comes to the 
equation describing cylindrical equilibria. We will refer to such a situation as an equilibrium 
collimation in the sense that at any z, the jet structure is the same as in an appropriate 
equilibrium cylindrical configuration. 

At \E' ^ \E', the plasma is significantly accelerated in the far zone so that the transfield 
equation is reduced to: 

^^2+^2^3^ ^2^2 ■ ^^^^ 

This equation describes the structure of the main body of the jet. One cannot give a simple 
physical interpretation of terms in this equation however, one can gain some physical insight 
considering regimes when different terms dominate. If the condition (15^ is fulfilled across 
the jet, one can neglect the derivative in z thus coming to a \Ef 3> ^E' limit of Eq. (HOl) . In 
this case, the jet as a whole is coUimated in the equilibrium regime. On the other hand, in 
some configurations (and anyway far enough from the axis) the condition opposite to (l39l) 
is fulfilled; then the term with the second derivative becomes dominant so that the equation 
is reduced to 

d^r d rfiii - 7)^ 

This equation could be dire ctly obtained assuming that the field is purely toroidal whereas 



the flow is purely poloidal (iLyubarsky fc Eichlerl I2OOII ). Then the flux freezing condition 
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dl]) yields B"^ — = {B^/'~^)'^. Substituting this relation into the transfield equation f|T^ 
and dropping the terms with Bp and v^, one comes, in the far zone, to Eq. P2l) . In 
this case, the function ^> could be considered as a specially normalized stream function, 
P7V oc (l/r)V\E' X 0, which describes the flow lines. The poloidal field. Bp, and the angular 
velocity, VL, become just auxiliary quantities formally defined by Eqs. ([5]) and ([6]). Note 
that Eq. (142|) does not change under the transformation \E' a^^f2; r] a^^rj; 

where a is an arbitrary number. As it follows from Eqs. (Q, and (fTOj) . neither poloidal 
velocity nor the electric field change under this transformation whereas Bp and Q could 
acquire any values. The flow with the poloidal field and azimuthal velocity neglected could 
be seen as composed from coaxial magnetic loops moving away and expanding together with 
the plasma. In this case, the difference between the hoop stress and the electric force is not 
counterbalanced by the pressure of the poloidal field. Taking into account that the electric 
field could not compensate the hoop stress completely (in the frame moving with the loop, 
the r component of the electric field is zero), one concludes that there is a residual force 
towards the axis of the flow. This does not mean that the flow immediately converges to the 
axis because in highly relativistic flows, the residual of the hoop stress and the electric force 
is small. In any case, we will refer to this situation as a non-equilibrium collimation. We will 
see that different regimes of collimation correspond to different acceleration regimes. 



4.3. Asymptotic Bernoulli equation and boundary conditions 

The transfield equation should be supplemented by the Bernoulli equation. We have 
already used this equation in the zeroth order in 1/r, Eq. (1271) . when simplified the transfield 
equation. However, one should be careful when using this equation in order to find 7 because 
7 turns out to be a small difference between two large terms if the flow is Poynting dominated, 
/i ^ 7. Therefore 7 could be found from the Bernoulli equation in the form of Eq. fl27j) 
only if a is not too large. Generally one should retain higher order terms and use, e.g., Eq. 
([2^1). Without loss of accuracy, this equation could be written as a cubic equation for 7 (e.g. 



Beskin et al.l fll998[ )) 



^ _ Q^r-'Bp n^r^Bp / _ 7^ - 7^ \ ^^^^ 

This equation is reduced to the zeroth order Bernoulh equation (127|) if one could neglect the 
expression in the right-hand side. This expression is small as compared with /i however, it 
could be neglected only when it is less than 7, i.e. only if 7^ ^ /i; (Qr)^ ^ /i . Note that 
7 = p} ^^ when the flow velocity is equal to the fast magnetosonic velocity (e.g., ICamenzind 



( 119861 )) so that one can find 7 from the zeroth order Bernoulli equation only beyond the fast 
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magnetosonic point. Well within this point, the Lorentz factor could be found from another 
limit of Eq. (03]): 

Generally one has to solve the cubic equation fj43l) so that there is no simple expression for 
7 valid in the whole far zone. This in fact means that the acceleration regimes inside and 
outside the fast magnetosonic point could be different. For example, in the split monopole 
wind, which represents a non-confined wind from a point source, the Lorentz factor grows 
linearly with the radi us until the fast ma gnetosonic point and then the acceleration becomes 



logarithmically slow (IBeskin et al.lll998l ). 



Note that the asymptotic transfield equation (138|1 is valid both outside and inside the 
fast magnetosonic point because it was derived only under the condition Qr ^ 1. In the next 
section, we show that when considering coUimated, Poynting dominated flows, one can avoid 
finding 7 from the Bernoulli equation. In this case, the acceleration regime does not change 
at the fast magnetosonic point therefore this point will not appear more in this paper. 

The asymptotic form of the boundary condition fl2U]) may be found by making use of 
Eq. (!29!) and taking into account that 7 ^ 7in in the outer parts of the Poynting dominated 
jet; this yields 

/ QrB„ \ ^ , , , , 

87rpext(^). (45) 



^ / >I'{r,2)=*o 

Taking into account Eq. (127|) . one could write this condition also as 



f2r7 



= 87rpext(2;). (46) 



5. The Poynting dominated flow in the far zone 

Let us first consider the structure of the Poynting dominated flow, /i ^ 7. Since the 
Poynting flux goes to zero at the axis, (see Eq. fl2^ ). this approximation is violated close 
enough to the axis, < \E'. Moreover, we will flnd that the flow is accelerated in such a 
way that the closer the fleld line to the axis, the earlier (at a smaller z) the flow kinetic 
energy approaches the total energy. Therefore a cr ~ 1 core is anyway presented within the 
Poynting dominated jet so that the results of this section could not be applied close enough 
to the axis. In Section 8, we flnd the structure of the flow close to the axis, which is smoothly 
matched, at a larger r, with the solution for the Poynting dominated flow. 
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5.1. The governing equation 

Here we study the structure of the flow at \E' ^ \[^, i.e., when the Poynting flux initially 
exceeded the plasma kinetic energy. In this case, we can use the asymptotic transfleld 
equation in the form fHTl) . 



As it was discussed in sect. 4.3, one cannot flnd a simple expression for 7 from the 
Bernoulli equation in order to substitute it into the transfleld equation and obtain a single 
equation for \E'. On the order hand, 7 could be easily found from the transfleld equation 
provided the shape of the magnetic surfaces, r(\E',z), is known. An important point is that 
in this case, an extra accuracy is generally not necessary because in the transfleld equation, 
7 is not presented as a difference of large terms. A special care should be taken only if the 
flow becomes nearly radial because the curvature of the flux surfaces is determined in this 
case by small deviations of the flow lines from the straight lines; this case will be specially 
addresses in sect. 7.2. In this and the next sections, we will neglect corrections of the order 
of 7//i to the shape of the flux line; then the Bernoulli equation ^17) is reduced to 

= r//i; (47) 

which could be considered, with account of Eq. flS^ . as an equation for r(\E', z): 

dv 

/xry— = rfil (48) 
The solution to this equation is presented as 



D{z)^{^y, $(^) = V2exp (^j^ 



(49) 



where D{z) is an arbitrary function. One sees that the structure of collimated, Poynting 
dominated jets is generally self-similar. Recall that this equation is valid only at ^ \E'; 
the solution will be continued to the axis in the Section 8. In any case, D is roughly the 
radius of the very inner part of the jet, ~ \E'. 

Close enough to the axis, one can use Eq. fl24l) for /i, which implies 

* = m 

This means the poloidal magnetic fleld becomes homogeneous, oc r^; dBp/dr = 0,well 
inside the jet, \I' <^ \I' <^ \E'o- Note that when flnding the expression fl24|) for fi, we assumed 
that the poloidal fleld is homogeneous near the axis so that this result is nothing more 
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than a consistency check. The same expression for is also obtained if one uses the 

general estimate fl26l) for /i. This is also because the coefficient 2 in fl26l) corresponds to the 
homogeneous poloidal field. Another coefficient would result in a power law function <I>(\E'). 
Such a strong dependence on /i arises because fi enters in the exponent. 

In order to find the function D{z), let us substitute Eq. poj) into the left-hand side of 
Eq. (HT!) and integrate the obtained equation between \l/ and \I'o: 



2D- 



dz 



2 /- 







^7/ *=vi> 



(51) 



Note that the region \E' ~ \E'o contributes mostly into the integrals therefore we could choose 
\l/ as the lower limit of integration even though the solution fHOl) is no longer valid there. 
One can also neglect the last term in the right hand side as compared with the first one 
because it could be checked a posteriori that the expression in the brackets grows with r. 
Making use of the boundary condition fH6|) . one reduces the right-hand side of this equation 



(52) 



to Svrr^pext = 47r$^D^Pext- Then one gets the equation for D{z) in the closed for 



d^D 



'to I 



*0 



firjd'^ 



dz^ i$ ^' y$ ^]2$2 

This equation could be written in the dimensionless form as 



-47r [$(^o)]V(^)^- 



where 



a 



Z = n{<ilo)z; Y{z) = a"^/^n(^o)^(^ 

4 fi'-f ' ; 



^^fir]d^j 



(3 = Anpo 



^^H'qd'i' 



(53) 

(54) 
(55) 

(56) 

Po=Pe.t{z = l/mo)); PoV{Z)=p{z). (57) 

This equation generahzes the equation for the jet radius obtained by Komissarov et al. ( 20081 ) 
as an order of magnitude estimate. We see that this equation is in fact asymptotically exact. 
Moreover, finding Y{Z) from this equation, one finds the full structure of the fiow. Therefore 
we will call Eq. fl53p the governing equation for Poynting dominated jets. 

Note that there is one to one correspondence between the terms in the governing equa- 
tion and in the original asymptotic transfield equation (j^Tl) . Namely the pressure term (the 
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last one) in Eq. fl53|l comes from the right-hand side of Eq. (HTl) whereas the first two terms 
correspond to the terms in the left-hand side. Following the discussion in sect. 4.2, one 
sees that the coUimation is in the equilibrium regime if the second term dominates the first 
one. Then one immediately finds Y{Z) = [f3V{Z)]~\ Of course neglecting the derivative in 
the equation, one looses solutions. The lost solutions just describe oscillations with respect 
to the equilibrium state. If the jet is not very narrow, the term becomes negligibly 
small; then the governing equation becomes linear. In this case the jet is collimated in the 
non-equilibrium regime. In Sect. 6. we present a more detailed analysis for the case of a 
power law profile of the external pressure. 



^The solution to the governing equation fl53l) could be presented ( jPolyanin &: Zaitsev 

20021 ) in the form Y = wy, where the auxiliary function, w, satisfies the linear equation 



Then the equation for y has the first integral 



(58) 



dZ 



(59) 



which could be immediately integrated once again. Now the general solution to Eq. (1531) is 
found as 

2" 



Y 



w 



1/2 



(60) 



In section 6, we present such a solution for the jet with a constant angular velocity confined 
by the external pressure decreasing as a power law. 



5.2. Finding the structure of the flow 

According to Eq. fH^ . the flux surfaces are self-similar in the Poynting dominated 
domain. Having found Y from the governing equation, one finds the shape of the magnetic 
surfaces as 

n(^o)r(^, Z) = a^/^$(^)y(Z). (61) 

Taking into account that the region \& ~ \E'o contributes mostly into the integrals in Eqs. 
f lSS]) and (jSni), one can estimate the coefficients in the governing equation as 
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In the last equality, we used the estimate (|26|1 . Substituting the obtained estimate for a into 
Eq. fl6Tl) . one finds 

r ~r(^o)^^(^o); (63) 

so that Y is of the order of the dimensionless outer radius of the jet. Making use of the 
estimate ^ ~ (l/2)r^i?p, one can write 



fi2 



(64) 



n(<I'o)r(<I'o)=l 



SO that j3 is of the order of the ratio of the external pressure to the magnetic pressure at the 
base of the fiow. 

In order to find the Lorentz factor of the fiow, one substitutes Eq. (1611) into the left 
hand side of Eq. fHTl) and performs integration between ^ and ^ to obtain 



(65) 



Retaining only the first term in the right-hand side, as it was done in Eq. (!52|) . one gets the 
relation for 7(\E', Z) in the closed form 



+ 



(66) 



Specifically for the periphery of the fiow, \I' = \I'o, one finds, with the aid of Eqs. fl53l) and 



7(^0,^) 



W 



^W{Z)Y{Z) 



-1/4 



(67) 



where 



With the aid of Eq. (137|) . one can write Eq. (166|) as 

r 



1 r 1 



AW 



-1 2 ^xlf 



5(vl/) 



$(^)r/(^)/i(^) 



-1 2 



(68) 

(69) 
(70) 



This e quation generalizes the equation obtained by lTchekhovskoy et al.l (120081 ) and by lKomissarov et al, 
(120081 ). Making use of the estimate (l26l) for /i, one finds that the coefficients A and B are 
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always of the order of unity. In the case of equihbrium coUimation, when the condition fl39l) is 
fulfilled, one ca n neglect the first term in the right-hand side, which yields the old-established 
( lBuckleylll977l ) acceleration regime 7 oc fir. In the opposite li mit of non-equilibri um coUi- 
mation, one comes to the scaling 7 oc ^/TZ/r found recently by lBeskin et al.l ( 2004 ). 



6. Poynting dominated jet with a constant angular velocity 

In this section, we apply the above general method to jets with the constant angular 
velocity, fi(^E') = const. In this case, one can conveniently use the dimensionless variables 

X = fir; Z = Vtz. (71) 

We also assume that the injection is homogeneous, ri{^) = const. 

Note that in this case, one can get simple relations assuming that the energy integral, 
/i, is described by the linear function fl24l) not only close to the axis but across the jet. Then 
the poloidal flux is homogeneous, see Eq. fl50|) . The coefficients a and /3 defined by Eqs. 
fl55|) and fl56|) . correspondingly, are reduced to: 



« = 3hTr ; (72) 









-A 




)■ 




^0 




6vrpo 






fi4^2 




-°o 



^ " fi4vi/2 - 52 ; (73) 

where Bq = fi^\I'o is the characteristic magnetic field at the light surface. Now the fiux 
surfaces are described by a simple formula 



X = 3^/^ i^—j Y{Z). (74) 
The coeffici ents in Eq. (1681) are r educe d toy4 = l/3;-B = lso that one gets the equation 



obtained by iTchekhovskoy et al.l (120081 ). In terms of Y and \E', the equation for the Lorentz 



factor ( 16611 could now be written in the simple form 

72 dz^ ^ y^' ^ ' 

One sees that close enough to the axis, the second term in the right-hand side dominates, 
which yields the acceleration regime 

7 = X (76) 
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If y ^ vZ, so that if collimation is not very strong, the first term could dominate in the 
main body of the jet, \1/ ~ \l/o. Then the Lorentz factor is determined by the curvature of 
the magnetic surface, namely 

7 = ^U^ =V^- (77) 



In any case, at \I' = \I'o one finds, with the aid of Eq. ( l53l) . 

31/4 

7(^0, Z) = , . (78) 

Below we present not only general formulae for the parameters of the flow but also simple 
estimates with the aid of Eqs. ([72]), ([HD and fITSD . 

Let the external pressure be decreasing as 

^=^- (79) 

Then the auxiliary equation (1581) is solved via the Bessel functions so that the general solution 
to the governing equation (153|) could be found analytically. Taking into account that the 
governing equation is valid only at large Z, one can use only the appropriate asymtotics of 
the solution. Since the asymptotics depends on the sign of k — 2, let as consider different 
cases separately. 



6.1. The case k < 2 
In this solution to Eq. (1551) is presented as 

w = ^J_L_ (y^Z'-^^' ] ■ (80) 
At a large Z, this function is reduced to 



w 



f_) cos5; S = ^^Z^-^I^-^^-. (81) 



TT \(3 I 2- K 2- kA 



Substituting this auxiliary function into Eq. (!60|) yields the general solution to the governing 
equation 



^ !2-KfZ^ 



TT \I3 



1/4 



— cos^ S + Ci I C2 cos S + sin S 

Ci V 2 — K 



1/2 

(82) 
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At Ci = (2 — K)/7r, C2 = 0, this solution is reduced to a power law 

1/4 



>^=(y) ; (83) 



which could be found directly from Eq. (1531) by neglecting the first term (iKomissarov et al. 



20081 ). One sees that with this solution, the first term in Eg. (15^ is much less, at k < 2, 
than the second one so that the coUimation occurs in the equilibrium regime. The general 
solution fl82l) also expands as Z'^^^ but very long wave oscillations are superimposed on this 
expansion, which means that the flow could oscillate around the equilibrium state. Such 
oscillations are possible if the jet was injected not in the equilibrium state. The amplitude 
of these oscillations could be found by matching to the near zone solution at Z ~ 1. The 
spatial period of these oscillations increases with the distance as Z'^/^. 

Taking into account that the governing equation becomes algebraic in the equilibrium 
regime, one can generalize Eq. ([HSD to the general pressure distribution provided the pressure 
decreases not faster than 2;"^: 

Y = (/5P)-'/' . (84) 

One can see that the jet expands while the confining pressure decreases. When the jet 
eventually enters the region with the constant pressure, the jet becomes cylindrical. 

The Lorentz factor of the flow is found from Eq. flBBl) . For the smooth expansion 
described by Eq. flS^ one can neglect the first term in the right-hand side, which yields 



One sees that in accord with the general analysis in section 5.2, the Lorentz factor of the flow 
is proportional to the cylindrical radius, 7 oc X. If the energy integral is a linear function 
of \E', Eq. fl2^ . which is anyway the case well within the jet, this expression is reduced 
just to 7 = X, see Eq. (1751) . This estimate remains valid also for non-power law pressure 
distributions when the jet shape is described by Eq. fl84l) . If the jet oscillates with respect 
to the equihbrium expansion, as is described by Eq. fl82l) . the Lorentz factor also oscillates 
with respect to that given by Eq. (ISBj) . 



6.2. The case k = 2 



In this solution to Eq. (l58l) is 



_ , ^/ZcosS; S = ^/3-l/A\nZ; /5 > 1/4; 
^-S z('+VT^)/2. (3<l/A. ^^^^ 
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Now the general solution to the governing equation is 

C0S2 S + Cl COS S + ^7=y| siu S 

1 + Cl {C2 - ^==^^^ 



Y 



^1 



,^1/2 



^(i/2)yT=4;a 



1/2 



1/2 



/? > 1/4; 



/3 < 1/4. 



^7) 



The /3 > 1/4 solution is similar to that for the k <2 case. At specially chosen con stants, 
C\ = (/5 — 1/4)^/^, C2 = 0, it is reduced to a pure power law (iKomissarov et al.ll2008l ). 

^1/2 

^ = (/5-l/4)V4' (88) 

whereas generally long wavelength oscillations are superimposed on the overall expansion. 
With the solution flHHl) . both terms in the left-hand side of the governing equation are com- 
parable so that this case is an intermediate between the equilibrium and non-equilibrium 
collimation. 



At /3 < 1/4, the solution is reduced, at large Z, to a power law (IKomissarov et al.ll2008l ) 

Y = CZ^ 



k 



1 



A(3)I2. 



9) 



Note that l/2<fc<lso that the flow is coUimated but slower than in the case j3 > 1/4. 
The constant C in this solution is not defined; it could be found only by matching to the 
near zone solution. If the flow was not coUimated at Z ~ 1, there should be C ~ 1. This 
solution could be obtained directly by neglecting the second term in the governing equation 
f l53|) . This means that at /3 < 1/4, the collimation is non-equilibrium. 

The Lorentz factor of the flow is found from Eq. fl66l) . At /9 > 1/4, one substitutes the 
solution fl55]) . which yields the relation 



21)2^2 



a 



(3-1/ A 



1/2 



" ^2.T //3-l/4^'/' 



a 



5 $2 



(90) 



which yields 7 oc v ^ oc X. In this relation, the terms in the right-hand side are comparable 
at \I' ~ \E'o. When decreases, the first term decreases faster therefore well inside the jet 
one can retain only the second term. Making use of Eqs. fISUl) and (jSSD, one finds 7 = X. 
So in this case the Lorentz factor of the flow is equal to the dimensionless cylindrical radius, 
which is the general property of the equilibrium collimation. 



When P < 1/4, one substitutes the solution (189!) into Eq. (166!) to yield the relation 



772^2 
2fiV 



^2(l-fc) 



(91) 
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At \& ~ \l/o, the first term in the right-hand side dominates the second one; this could be 
easily seen from Eq. fl75|) , which is the approximate form of Eq. fl66l) . Then one finds 

-1/2 



( 2/?v^ j 7]^i^^dm\ Z^~^- (92) 

which reproduces, in this specific case, the scaling 7 oc \jTZ/r common to the non-equilibrium 
collimation. For ^ given by the lineat function fl2^ . this relation is reduced, with the aid of 
Eqs. ([50D and ([72]) to 

Note that in this case, the Lorentz factor increases, at a fixed Z, towards the axis. The 
Lorentz factor increases until at small enough \Ef, the second term in the right-hand side of 
Eq. (jnH) becomes dominant, which means that close enough to the axis, the jet is in pressure 
equilibrium. In this region, one has 




^ = a''''C\^Z' = X- (94) 



as in any equilibrium flow. So at any fixed Z, the Lorentz factor increases outwards from 
the axis while the flow is in the pressure equilibrium and then decreases outwards. 

The transition from the non-equilibrium to the equilibrium zone occurs at 

1 

^ " T^c^z^' ^^^^ 

when the two terms in the right-hand side of Eq. fl^T]) become equal. Transforming to the 
coordinate space with the aid of Eqs. (161 p and (1891) . one sees that the transition occurs at 
the line 

The Lorentz factor of the flow increases as 7 = X while the flow remains within the line 
( !96|) whereas after the flow enters the non-equilibrium zone, the acceleration proceeds slower, 
according to Eq. 



6.3. The case k>2 

In this solution to Eq. (l58l) is presented as 
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The large Z asymptotics of this function corresponds to the small argument limit of the 
Bessel function, 



J Ax) 



X 



r(z/ + 1) V2, 

where F is the gamma-function. Then the auxiliary function, w, goes to a constant 

1/(k-2) 



1 



W 



VP 



whereas the general solution (1^ goes, at large Z, to a linear function 



(98) 



(99) 



Y 



K-2 



K-2 



VP 



1/(k-2) 



Z; 



(100) 



so that the flow becomes radial at large distances. 



One sees from Eq. fllOOp that the flow could be coUimated, Y <^ Z, if k, only slightly 
exceeds 2 or/and P is large. Note that in this case, the small argument limit of the Bessel 
function, x -C 2, which yields a linear asymptotics for Y, is achieved only at a very large Z. 
For example, if k = 2.5, the above limit is achieved only at Z ^ 16/?^. In order to see what 
happens at a smaller Z, let us assume that 



VP 

K-2 



> 1. 



(101) 



Then the argument of the Bessel function in Eq. (1971) is large at 

2^n2/(-2) 



Z < 



K-2 



(102) 



In this case, one can use the large argument asymptotics of the Bessel function in Eq. (|9 
which leads to (cf. Eq. (|82D ) 



w 



TT 



Y 



K 



TT 



2 /Z^ 

P 



1/4 



COS S: 



S 



^VP ,1- 



1 



— COS S + Ci { C2 COS S 



IT 



K IT 

K -24'' 



■ sin S 



K 



1/2 



(103) 



(104) 



So the solution obtained for the case k < 2 could be extended to k slightly above 2 but only 
in a limited range of Z. As in the case k, < 2, the solution describes smooth expansion 



Y 



2K \ 1/4 

T 



(105) 



- 25 - 



only if the constants are specially chosen, 



Ci = - 2)/7r; C2 = 0. 



(106) 



Generally long wavelength oscillations are superimposed on the overall expansion. 

One sees that if the condition (1 10 II) is fulfilled, the flow is collimated according to Eqs. 
(11041) or (11051) in the region (11021) : at a larger Z, the flow becomes radial preserving the 
acquired collimation angle. In Fig. 1, solutions to the governing equation are shown for 
K = 2.5, (3 = 5. These solutions are obtained by numerical integration of Eq. ( [53|) . which is 
easier than numerical evaluation of integrals in Eq. (M?]) . The smoothly expanded solution is 
shown by solid line whereas dashed line represents a solution with oscillations in the region 
(11021) . All the solutions go to a linear function at large Z where the flow is already well 
collimated. The final collimation angle could be obtained from Eq. (llOOp as B = Y/Z. 
Choosing the constant Ci from Eq. (11061) corresponding to the smoothly expanded solution 
(11051) . one finds the final collimation angle as 



One sees that the collimation angle rapidly increases with increasing k and decreasing (3; 
= 0.01/(3^-^ atK = 2.2, e = 0.2//? at k = 2.5 and = 0.56/v^ at k = 3. This means that 
if the external pressure decreases faster than V oc Z^^, a narrow jet could not be produced 
unless the flow has already been collimated in the near zone, Z ~ 1. 

Note that the solution (11051) corresponds to the equilibrium collimation because it could 
be obtained by neglecting the first term in the governing equation. In order to figure out the 
collimation type of the radial flow (llOOp . one has to find the curvature of the field surface. 
Expanding the solution to the governing equation (153|) to higher order terms in 1/Z (this 
could be easier done by making expansion in the equation than by expanding the general 
solution (160|) and (1971) ). one gets 



where ci and C2 are constants, which could be expressed via the constants Ci and C2 in the 
general solution (1^ . Note that Ci is in fact the final opening angle of the jet, Ci = 0; below 
we will use the expression (11071) for Ci. One sees from Eq. (llOSp that at k > 4, the curvature 
is independent of the parameters of the external pressure, (3 and k, which means that the 
flow is not confined by the medium with such a sharply decreasing pressure. Below we do 
not consider this case. At k < 4, the ratio of the first to the second term in the governing 




(107) 




2< K<A; 
K> 4; 



(108) 
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Fig. 1. — The shape of the jet, Y{Z), for k — 2.5; (3 — 5. The solution without oscillations 
is shown by solid line; dashed line shows a solution with oscillations. Dotted lines show 
asymptotics, X oc Z**/^ = Z^^^ and X oc Z , correspondingly. 
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Fig. 2. — Lorentz factor of the flow shown by sohd hnc in Fig. 1 (expansion without 
oscillations). Solid line shows the Lorentz factor of the flow at the boundary of the jet, 
\E' = \I'o; the dashed line is for the Lorentz factor at the flux surface \1' = 0.2\t'o. Dotted lines 
show asymptotics, 7 oc Z^l'^ = Z^l^ and 7 oc Z'^'^~'^^l'^ = Z^^*, correspondingly. 
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equation fl53l) is 



- = pe'Z'--; (109) 



which means that the flow is in the non-equihbrium regime at 

Z>Z,= {PQ')-'^^'-^^ (110) 

For K = 2.5, the transition occurs at Zi = 70/3^. 

The Lorentz factor of the flow is obtained, as before, from Eq. fl66l) or, for the hnear 
^(^), from Eq. (I75|l . At the condition (11021) . when the flow is colhmated according to 
Eq. (11051) . the second term in the right-hand side dominates the first one (the equilibrium 
colhmation) so that one comes again to Eq. fl55]) . which is reduced, at \l/ <C ^'oj just toj = X. 
At Z > Zi, where the flow is aheady nearly radial, the first term becomes dominating, which 
yields 



7 



2/3 J r]f^^^d^\ (Ill) 



This is in accord with the general non-equilibrium scaling 7 oc >/TZ/r. For linear dependence 
of on \E', this relation is reduced, with the aid of Eqs. (I50l) and (I72|) . to 



7 



3^ /^^(«-2)/2_ (112) 



In Fig. 2, we show the Lorentz factor of the smoothly expanded (without oscillations) flow at 
the same parameters as in Fig. 1. One sees the transition from a relatively rapid acceleration 
in the equilibrium regime, 7 oc Z'^^'^ = Z'^l^ to the slow non-equilibrium acceleration, 7 oc 

An important point is that in the non-equilibrium regime, the Lorentz factor increases 
with decreasing \1/ so that the flow is faster inside the jet than at the periphery. At small 
enough the second term in Eq. ( 1751) could become dominating, the transition occurring 

at 

\ z) 

Recall that the non-equilibrium zone appears only bX, Z > Z\. In the coordinate space, the 
transition occurs at 

1/4 



(113) 



Inside the zone bounded by this surface, the flow is accelerated as in the equilibrium case, 
7 = X. 
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One has to stress that at Z > Zi, the Lorentz factor of the flow is determined by 
the curvature of the flux surfaces (non-equihbrium regime) even though the flux surfaces 
are nearly conical so that the Lorentz factor depends on small deviations from the conical 
shape. In this case, accuracy of the governing equation could become insufficient in order to 
find the Lorentz factor of the fiow. The governing equation was obtained by neglecting the 
kinetic energy term in the Bernoulli equation fl271) therefore the shape of the fiux surfaces, 
r(^, z), is determined to within a factor of 7///. In the case of the equilibrium coUimation, 
7 oc r, this accuracy is sufficient while the jet remains Poytning dominated. In the case of 
non-equilibrium coUimation, the Lorentz factor goes as \/TZ/r and if the flow lines become 
nearly straight, the curvature could be determined by the neglected terms of the order of 
7//i. In this case, the scalings fillip and (11121) cease to be valid when the flow is still Poynting 
dominated. We address this issue in sect. 7.2. 



6.4. Comparison to previous works 



Recently magnetic accelerat ion of externally con f ined j ets w as carefully studied, bot h 
numerically and ana l ytical ly, by iTchekhovskoy et al.l (120081 ) and iKomissarov et al.l (I20081). 



chekhovskoy et al.l (120081 ) used the force-free approximation whereas 



Komissarov et al. 



(2008) solved the full set of relativistic MHD equations. For analytical estimates. Ixchekhovskoy et al, 
(120081 ) assumed that the shape of the flux surface is a power law and then found from the 
transfleld equation the appropriate Lorentz factor ar id the external pre s sure. They obtained 
the scaling (183!) and claimed that it is un iversal. IKomissarov et al.l (120081 ) based on the 



asymptotic transfleld equation obtained by Vlahakia i 2004 ) . which is equivalent to our Eq. 
(IMl) . Order of magnitude estimate of terms in this equation led them to Eq. ( l53l) for the 
jet radius. Analyzing this equation, they revealed that the scaling (15^ is valid only for 
K < 2. Foi K = 2, they obtained the scalings ([88]) and ([89]) for /5 > 1/4 and /5 < 1/4, 
correspondingly. For k > 2, they obtained the radial asymptotics. 

Our approach generalizes these flndings permitting the asymptotically exact solutions 
describing the full structure of the get. Going beyond the simplest power law scalings also 
permit us to flnd important new qualitative features of the flow. In particular, we found that 
in the case of the equilibrium coUimation, k < 2, oscillations could be superimposed on the 
general expansion of the jet. For k > 2, we see the transition between the equilibrium and 
non-equilibrium regimes, which could not be described by a power law scaling. Namely, if k 
only slightly exceeds 2, the flow is coUimated according to the equilibrium law ([83]) but only 
till some limiting distance beyond which the flow becomes radial preserving the acquired 
coUimation angle. The larger [k — 2), the earlier (at a smaller distance) the flow becomes 
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radial so that at k > 3 the flow is practically radial from the very origin. 



Tchekhovskoy et al.l (120081 ) numerically simulated jets with different profiles of the ex- 
ternal pressure. They found excellent agreement with the scaling fl5^ aX k = 2. For k = 2.5 



they reported noticeable deviations from this scaling at large distances; we suppose that they 
observed the transition to the radial flow. At last for k, = 2.8, they found a wide conical jet. 
These results agree with our conclusions. 



7. Terminal Lorentz factor and collimation angle. 

It was shown in the previous section that if the Poynting dominated outflow is con- 
fined by the external pressure, the flow is coUimated and the plasma is accelerated so that 
eventually the kinetic energy of the plasma could not be neglected any more. An important 
point is that the closer to the axis, the earlier (at a smaller z) this happens. Note that at 
\1/ < the Poynting flux is relatively not large from the very beginning. In this section, we 
address saturation of the acceleration and estimate the terminal Lorentz factor and collima- 
tion angle. For the estimates, we assume that the energy integral is described by the linear 
function fl2^ : then the parameters of the flow are given by Eqs. fl50l) . fl72|) and (I73|l . It is 
also convenient to introduce the maximal achievable Lorentz factor, 

7max = /^(^o) ~ 7m^o/^; (115) 



which is just Michel's magnetization parameter (jMichellll969l ). 



The results of the previous section have been obtained in the limit 7 <^ which means 
that the shape of the flux surfaces has been found with the accuracy of 7//X. We showed 
that the flow is accelerated as 7 ~ X in the case of equilibrium collimation, i.e. under the 
condition fl39p . and as 7 oc ^/TZ/r in the opposite case. A small error in the shape of the 
flux surface yields the error of the same order in the equilibrium scaling 7 oc X therefore 
this scaling could be safely extrapolated up to 7 ~ /i ~ ji^'^/'^q, i.e. until the flow ceases 
to be Poynting dominated. In the case of non-equilibrium collimation, one has to analyze 
how corrections of the order of 7//X could alter the curvature, 1/TZ = ~(Pr/dz^. Small 
corrections to the shape of the flux surface could signiflcantly modify the curvature if the 
flux surfaces are close to cones, i.e. if r{Z) is close to a linear function. This indeed happens 
when the conflning pressure decreases faster than z'"^ . Therefore we consider separately the 
cases K < 2 and k > 2. 
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7.1. The case k < 2; transition to a ~ 1. 

In the case k < 2 the Lorentz factor of the flow increases according to the equilibrium 
law 7 = X. Corrections of the order of 7/yU <^ 1 could not affect significantly this scaling 
therefore the flow is accelerated up to 7 ~ yU ~ 7max^/^o; this occurs at X ~ 7max^/^o- 
Making use of Eq. (1711) and fl55]) . one finds the corresponding distance as 



(- 

3 V^o 



1/k 

(116) 



Reverting this expression, one finds the boundary of the moderately magnetized, a ~ 1, core 
as 

^core /3 ^^^^^ 



^0 V /3 7, 



y2 • 

/max 



In the coordinate space, the boundary of the moderately magnetized core is found as 

Xrnrp = \ ^ • (118) 




One sees that the closer to the axis, the smaller the distance where the acceleration 
saturates so that a mildly magnetized, a ~ 1, core expands with the distance occupying a 
progressively larger fraction of the jet body until the a ~ 1 state is achieved across the whole 
jet. This happens at the distance 



3 ""^ 

from the origin. At this distance, the collimation angle, G = dY/dZ, is 



(119) 



q(4-k)/4k 

e = - — -^^-(4-)/-. (120) 

'max V-^^"/ 

When the flow ceases to be Poynting dominated, the collimation angle decreases further. It 
will be shown elsewhere (Lyubarsky, in preparation) that if k < 2, the flux surfaces become 
cylindrical at the infinity and the Poynting flux is totally transferred to the kinetic energy. 

If the density decreases not as a power law (but slower than z""^), one can write general 
estimates making use of Eq. (IMl) and flTHj) . Specifically one finds that the pressure should 
decrease at least by a factor of /?7max order for the Poynting flux to be converted into the 
kinetic energy of the flow. If the environment at large distances from the compact object 
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has the finite pressure pamb > Po/(/57max) ~ -^o/(67i"7max)? the fiow becomes cyhndrical and 
the Lorentz factor is saturated at the value 

/ ^2 \ 1/4 

The above estimates are illustrated by a sketch in Fig. 3. The curve 1 shows the 
distribution of the Lorentz factor across the jet at a not very large distance from the origin 
where the fiow is still Poynting dominated everywhere with except of the region \E' < \E'. 
Further out of the origin, the a ~ 1 core expands within the jet. The curve 2 shows the 
Lorentz factor at some intermediate distance where the core is already developed but the 
main body of the jet is still Poynting dominated. The curve 3 sketches the distribution of the 
Lorentz factor at the distance Z ^ Z^j when the whole jet ceases to be Poynting dominated. 

The above estimates could also be applied to the case k = 2, /?>l/4 when the fiow 
follows scalings obtained for k < 2; one has just substitute k by 2 and /? by /? — 1/4. For 
example, the boundary of the moderately magnetized core is now presented as 



whereas the distance where the jet ceases to be Poynting dominated is 

-tLx- (123) 



On the contrary, if k = 2 and /5 < 1/4, the fiow exhibits quite different behavior 
because it is coUimated in the non-equilibrium regime. In such a flow, the Lorentz factor 
grows oc X only if and while the flow line remains close enough to the axis where the flow 
is in the equilibrium. After the flow line crosses the boundary between the equilibrium 
and non-equilibrium zone given by Eq. ( 196|) . the acceleration proceeds slower, the Lorentz 
factor being given by Eq. ( 193|) . In the Poynting dominated domain, the Lorentz factor is not 
monotonic across the jet; at a flxed Z, it increases with the radius within the equilibrium zone 
and decreases outwards in the non-equilibrium zone. Therefore at any Z, the Lorentz factor 
reaches the maximum value somewhere within the jet. Comparing the Lorentz factor with 
/i(\E'), one flnds that at a distance Z from the origin, the flow remains Poynting dominated 
only outside the boundary 

^ ^ r ^?.C-TZ.Z-^ Z < [7n.ax/(3V^/5V4C2)] '/^''-'^ ; 

^0 \ 3V6 (/?C^Lx)"'^' Z > [7max/(3l/^/5V4C2)] '^^''-'^ . ^ ' 
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Fig. 3. — Distribution of the Lorentz factor across the jet in the case k < 2; not to scale. 
Thick hne shows the distribution of the total energy, /i(\E'). Each thin line shows the distri- 
bution of the Lorentz factor at some distance from the origin; they are labeled according to 
the distance, i.e. the curve 1 is the closest to the origin, the curve 3 is the farthest. 
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or, in the coordinates, 



One sees that the acceleration in the central part of the jet should be rapidly saturated 
because the Lorentz factor achieves the maximal value 7 ~ Slow acceleration in 

the outer, non-equilibrium part of the jet continues until the Lorentz factor approaches 
The closer the flow line to the axis, the earlier (at a smaller Z) this happens so 
that in this case also a moderately magnetized core expands within the body of the jet. At 
Z > [7max/(3^''^/3^/^C^)] ^'''■^'^ ^\ a maximal Lorentz factor is achieved at the boundary of 
the core because in the non-equilibrium jet, the acceleration is faster at the flow lines closer 
to the axis. The whole jet ceases to be Poynting dominated at the distance 

^conv= 3-l/^V^C7max • (126) 



The corresponding coUimation angle is 



e = 4-^- (127) 

VP7max 



The distribution of the Lorentz factor across the jet in the case k = 2; [3 < 1/4 is sketched 
in Fig. 4. The curve 1 shows the Lorentz factor not far from the origin. An internal part of 
the jet is collimated in the equilibrium regime and the Lorentz factor increases outwards from 
the axis. In the main body of the jet, the collimation is non-equilibrium and the Lorentz 
factor decreases outwards so that the Lorentz factor is a maximum at the boundary between 
the equilibrium and non-equilibrium zones. As the jet propagates, the flow in the internal 
parts reaches the a ~ 1 state and stops accelerating. The curve 2 shows the distribution 
of the Lorentz factor at the distance Z > [7max/(3"'^/'^/9^/^C^)] ^''^^'^ where the flow outside 
the core is non-equilibrium so that the maximal Lorentz factor is achieved at the boundary 
of the core. The curve 3 shows the distribution of the Lorentz factor at Z ~ ^conv when the 
whole jet ceases to be Poynting dominated. 



7.2. The case k > 2; transition to logarithmic acceleration. 



At K > 2, the flow is collimated only if the condition (110 II) is fulfilled. In this case, the 
flow becomes conical still being Poynting dominated; the final collimation angle is given by 
Eq. (11071) . In the conical part of the jet, the Lorentz factor grows according to Eq. (I112p . 
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Fig. 4. — The same as in Fig. 3 but in the case k, = 2, f3 < 1/4. 



- 36 - 



which corresponds to the general non-equihbrium scahng 7 oc \jTZ/r. An important point is 
that the curvature of the flow hues is determined by small deviations from the straight line, 
see Eq. (11081) . therefore Eq. (lll2p is valid only if the shape of the flow lines could be found 
from the governing equation with the necessary accuracy. Let us consider more carefully the 
jet with nearly straight flux surfaces. 

In the previous sections, we found the jet structure neglecting 7 as compared with /i in 
the Bernoulli equation (1271) . Then the shape of the flux surfaces is presented as (see Eqs. (pTl) 
and (CID) X(*, Z) = a^/^^{<i/)Y{Z) = 3^/^^/¥/W^Y{Z), where Y{Z) satisfles the governing 
equation. We can flnd limits of applicability of Eq. (I112p for the Lorentz factor of the flow 
by substituting this equation into the Bernoulli equation fl27j) and flnding the corresponding 
corrections to the shape of the flux surfaces. Eq. (11121) is valid while the curvature due to 
this corrections remains small as compared with the curvature obtained from the solution 
(11081) of the governing equation. 

Let us present the shape of the flux surfaces as (cf. Eq. (!6T!) ) 

X = a^/^$y(l + (5); (128) 

where 5(\1', Z) <^ 1 describes corrections to the shape of the flux surfaces due to a non-zero 
■j/fi. Substituting this into (1771) and linearizing with respect to small 5 and 7///, one gets 

^ - ^ (129) 

Assuming for simplicity that the energy integral is described by the linear function flM|) . one 
writes in the dimensionless form 

d6 7 



where 



With 7 from Eq. (I112p . one flnds 



^ = ^/^o- (131) 



3l/4v(K-2)/2 

(132) 



Substituting this expression into Eq. (11280 and differentiating twice with respect to Z, one 
finds the curvature of the flux surface as 

^ ^ ,1/4 [Yffl_ 3V^/.(/. - 1) \ 

dZ^ y^0\dZ^ 4v^7max^3/2^2-''/V ■ ^ ' 
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Here we take into account that Y ^ QZ. The Lorentz factor of the flow could be determined 
from the solution to the governing equation only if the second term in brackets is small as 
compared with the first one. Finding d'^Y/dZ'^ from Eq. fllOSp . one sees that this is the case 
only at distances smaller than 



47max 



2/[3(k-2)] /^yNl/(K-2) 



_3V4fi;(K- 1)_ 

At this distance, the flow acquires the Lorentz factor 



(134) 



Note that this Lorentz factor is the same for all flux surfaces whereas Zt increases towards 
the periphery of the jet. This is because in the non-equilibrium regime, the acceleration rate 
decreases outwards from the axis. 

Let us now find the Lorentz factor of the flow at Z > Zt. With this purpose, one has 
to solve the transfield and the Bernoulli equations without neglecting 7 in the Bernoulli 
equation. In the case of interest, the coUimation is non-equilibrium therefore one can take 
the transfield equation in the form (H2il . As the flow lines are nearly straight, we can look 
for the solution in the form (11281) with Y{Z) = QZ. We again assume for simplicity that 
the energy integral is a linear function flM|) : then $ and a are given by Eqs. fl50l) and fl72|) . 
correspondingly. Now the transfield equation is written in the limit 5<^1, 7/yU<^las 

_y3e'z(2- + ^)=---. (136) 

Linearization of the Bernoulli equation in small 6 and j/fi yields Eq. (11301) . 

Eliminating 6 from these two equations (by differentiating Eq. (I136P with respect to S 
and substituting Eq. (I130p ). one gets a single equation for 7: 

('^az + ^ az? j = 1 + '7as + ^ ^ as (yfas j ■ 

As an initial condition, one can take 7 = 74 at Z = (the above estimates give in fact 
7 ~ 7f at Z ~ Zt). Note that 7^ is independent of S. Assuming that beyond Z^, the solution 
is also independent of S", one finds with the logarithmic accuracy, i.e. in the limit InZ ^ 1, 

7= ^^^^|^(lnCZ)l . (138) 
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One sees that with the constant C = 1/Zt, this function goes to 7 ~ 7t at Z ~ and still 
satisfies Eq. (11371) with the logarithmic accuracy (because it only logarithmically depends 
on S, via Zf). So the final solution at Z ^ Zt (in fact at In Z/Zf ^ 1) is written as 




7=r^^ln^ • (139) 



According to this solution, the fiow in fact stops accelerating beyond the distance Zt so that 
one can use 74 as an est imate for the terminal Lorentz factor . This conclusion matches with 



the well-known result (ITomimatsul 1 19941 : iBeskin et al.l Il998l ) that the radial, non-confined 



wind is accelerated only till 7 ~ 7max and then the Lorentz factor grows only as (Ini?)^/^. 



One has to stress that according to the boundary condition (I46il . the flow at the bound- 
ary is accelerated till 7max provided the external pressure falls to zero. Therefore close enough 
to the boundary of the flow, the Poynting flux is efficiently converted into the kinetic energy. 
Specifically for the power-law pressure profile ( j79|) . the boundary condition ( l46l) yields in the 
limit 7 < 7max (El 

^i'^o,Z) = ^Z^^-'y'; (140) 

which recovers the scaling (I112p . On sees that at the boundary of the fiow, the scaling (11121) 
remains valid until 7 ~ 7max even though in the main body of the jet, the acceleration is 
saturated at 7 ~ 7^. In order to find the width of the boundary region where the acceleration 
proceeds beyond ■jt, let us substitute Eq. (11401) into Eq. (I137P and estimate d'^'j/dS'^ ~ 
■-y/{ASy necessary to satisfy the equation. This yields 

AS.^.f£i(M)~. (141) 
\E'o \ Z J 

One sees that after the fiow reaches 7 ~ 7^ at Z ~ Z^, the acceleration proceeds further only 
in a narrow region close to the boundary. 

The results of this subsection are illustrated in Fig. 5. The curve 1 shows the distribution 
of the Lorentz factor at some distance Z < Z\ where the jet is coUimated in the equilibrium 
regime. The curve 2 corresponds to a distance Z > Z\ where the main body of the jet is 
in the non-equilibrium regime so that 7 has a maximum inside the jet, at the boundary 
between the equilibrium and the non-equilibrium zones. The curve 3 shows the distribution 
of the Lorentz factor oX Z = Zt{^^ for some \E' < \l'o. At this distance, the acceleration is 
saturated in the internal part of the jet. At last the curve 4 shows the Lorentz factor at 
Z > Zt(\E'o) where the acceleration is saturated in the main body of the jet and proceeds 
further only in a narrow boundary layer. 
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Fig. 5. — The same as in Fig. 3 but in the case k> 2. 
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Note also that the above estimates assume that the flow remains Poynting dominated, 
i.e., that 7^ < 7max^/^o- In the opposite case, the acceleration is saturated at 7 ~ 7max^/^o- 
In particular, if 7^ > 7max; the total Poynting flux is efficiently converted into the kinetic 
energy. The last condition could also be written as 7max0 < 1? which means that the flow 
remains causally connected in the sense that a signal sent in the transverse direction (in 
the proper plasma frame) could cross the jet for the proper time zj'^. Recall that the same 
condition is satisfied when k <2 therefore in any case, the whole flow is accelerated till a ~ 1 
only if it remains causally connected. If the acceleration is saturated at 7 ~ 7^ < 7max, the 
flow is causally disconnected, 7t9 ~ (7max0)^''^ > 1- By this reason, in particular, the flow 
do not "feel" the boundary any more so that the acceleration stops everywhere except of a 
narrow boundary region. The linkage betwee n acceleration and causal connection of the flow 
was also noted by iTchekhovskoy et al.l (120091 ). The loss of causal connection also implies the 
global MHD stability of such jets because global instabilities (e.g., the kink instability) has 
no time to develop. 

Note also that 6 is determined only by the external pressure profile whereas 7max is 
determined by the parameters of the outflow so that these two quantities are independent 
and any relation between them is possible. However, 7^6 could hardly ever be very large 
in real systems because this quantity depends on the parameters in the power 1/3. On 
the other hand, such a flow could be accelerated further when and i f the Poynting flux is 



2002 




dissipated (jThompson 



Kirk fc Skj38raasenll2003l ) 



One should also note in this connection that according to a widely accepted view, the 
observed achromatic breaks in the GRB afterglow light curves occur when G7 becomes 
approximately unity. Since the afterglow is attributed to the decelerating jet, this implies 
that 7G was larger than unity in the prompt phase. One sees that the required property 
could be achieved in the MHD scenario if the confining pressure decreases with the distance 
something faster than 



8. The core of the jet 

It was shown in the previous section that a moderately magnetized core occurs near 
the axis of a Poynting dominated flow so that the solutions obtained in the section 6 could 
not be continued to the axis. In this section, we find the structure of such a core smoothly 
matched with the structure of the Poynting dominated flow in the main body of the jet. 

The flow near the axis is described by Eq. fHUj) . which should be complemented by the 
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Bernoulli equation in the form ([2 

r)X 

^(/i-7)^ = f^^X (143) 

This set of the first order differential equations for X{^) and 7(^1^) should be solved at the 
condition X(0) = and matched, at large \E', to the solution in the main body of the flow. 
For example, if the main body of the flow is Poynting dominated, the solution to Eqs. (11421) 
and (I143P should be matched with the solutions obtained in Section 6. The dependence on 
Z enters only via this matching. 

Close to the axis, the energy integral has the form of a linear function (I24p . Introducing 
the variables 

s = i + ^- e = — ; r = -L; (144) 

^ 7in 7in 

one reduces Eqs. (I142p and (11431) to dimensionless form 

Yds 2(i + e)(s-r)' ^ ' 

^ = ^ . (146) 

ds 2(s-r) ^ ' 

Before presenting the numerical solution to these equations, let us investigate them qualita- 
tively. 

Near the axis, ^ <^ 1, the solution is 

s = l + Ce-, r = 1 + ^Ce-, (147) 

where C is a constant. This means, in particular, that the poloidal magnetic field is homo- 
geneous at \E' <^ \E'. 

Far from the axis, s ^ 1, ^ 1, the flow is accelerated, F ^ 1, so that Eq. (11451) is 
reduced to 

sdT sF^ ^ ^ 

r d7 = ^ - 

The set of equations (11461) and (11481) is invariant with respect to the transformation s As; 
,^ ^ A^; F ^ AF so that the equations could be integrated. Namely, introducing the variables 

M = ^; a=^-l; (149) 
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(note that a thus defined is indeed the ratio of the Poynting to the kinetic energy fiuxes) 
yields the set of equations 



2(TS^ = (1 + (t)V; 
as 

_as du 



u ds 



a 



which has the first integral 



(a -I] 



Ci. 



The general solution is written as 

s = C2 [I - a + V -CiJ [I - (J - V -ci 

s = C2 [ci + ((T - 1)^] exp 



2 o-l 
arctan 



ci < 0; 
ci > 0. 



(150) 
(151) 

(152) 

(153) 
(154) 



/Ci ^Ci 

The solutions with ci < describes the fiow with a growing with the radius, and therefore 



with until it reaches a constant o"o = 1 — y— Ci < 1. Therefore this solution represents 
the structure of the core in low-cr jets. Transition to cr ^ ctq is described by the expression 



a = ao-[2(l-ao)]( 



(2-ao)M, /C2\-(i-o)/-o 



(155) 



Making use of Eqs. fll49p and (11521) . one writes this asymptotics in the original variables as 



1_ 

7in 



1 + O-Q ^ 



X , , . _ 
— = [2(1 - ^o)] 

7in 



-1/cto (1-cto)/2o-o 



(l+cro)/2<TO 



(156) 



The solutions with ci > become Poynting dominated far enough from the axis. In the 
limit cr ^ 1, one finds 



M = 1; s = C2<J exp 



TT 



or, returning to the original variables. 



7 = X; ^ = C2 exp — — 

* V V^J V7in 



xy 



(157) 



(158) 



One sees that far from the axis, the poloidal field becomes homogeneous and the solution 
is smoothly matched with the solution for the Poynting dominated domain, Eqs. (ISUI) . (IMl) 
and (175]). 



When ci is not small, Eq. (11571) implies that a is large at large s so that the fiow is 
Poynting dominated everywhere except of the region < \E'. When ci is small, the Poynting 
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dominated domain arises only very far from the axis, at s ^ exp(7r/^/ci) ^1. In the 
intermediate region, 1 <^ s <^ exp(7r/y^), the solution (I154p is reduced to an intermediate 
asymptotics 

a = l + -ln-. (159) 
2 ci 

This means that in the core of the Poynting dominated jet, a is close but remains larger 
than unity. It could become less than unity only when the whole jet ceases to be Poyniting 
dominated. In the original variables, the solution in the intermediate zone is 

~ = ; 7 = i^c/ X. (160) 

^ 7in ^ 

Numerical solutions to Eqs. (11451) and (11461) are presented in Fig. 6. For all asymptotics 
to be seen clearly, we plotted the curves in a very large scale. The curves are labeled by 
an appropriate constant C in the left boundary condition (I147p . When C > 1.5, the flow is 
Poynting dominated everywhere except of the region \E' < \E'. The poloidal magnetic field 
is homogeneous as it should be in the Poynting dominated flow with the energy integral 
/i(\E') given by Eq. (^^. In the case 0.38 < C < 1.5, the solution goes to the Poynting 
dominated asymptotics (I158P only at large enough distances from the axis; Bp goes to a 
constant in this zone. Between \[f ~ \[' and the Poynting dominated zone, the solution is 
roughly described by an intermediate asymptotics (11591) and (I160p . The poloidal magnetic 
field varies roughly as oc 1/X in this zone, which means that the toroidal field, = XBp, 
and the Poynting flux remain roughly constant. At C < 0.38, the solution is described, at 
^ ^ 1, by the asymptotics (11560 . corresponding to cr = const < 1. In these solutions, the 
poloidal magnetic field decreases faster than 1/X. These solutions describe the core of the 
jet at the stage when most of the Poynting flux is already converted into the kinetic energy. 

Any solution to Eqs. (I145P and (I146P with the initial condition (I147P describes the 
transverse structure of the jet at some Z. Generally a decreases with Z therefore the curves 
in Fig. 6 describe the Z development of the jet "upside down", i.e. the upper curves describe 
the transverse structure of the jet at smaller Z. The moderately magnetized, a ~ 1, zone 
occupies initially only the region ~ As the distance grows, the a ~ 1 zone extends to 
a larger 

The structure of the core of the Poynting dominated jet may be found in any specific 
case as follows. First one finds the shape of the flux surfaces in the Poynting dominated 
region, Y{Z), as it was described in the previous section. Then the transverse structure 
of the core at any Z is described by a solution to Eqs. (11450 and (I146P satisfying the left 
boundary condition (I147p and matched, at large s, with the inner solution for the Poynting 
dominated flow, Eqs. fl50l) . fl6T|) and fl72|) . The matching is reduced to finding an appropriate 



-44- 





Fig. 6. — The transverse structure of the core of the jet. Shown are the ratio of the Poynting 
to the kinetic energy flux (a), distribution of the poloidal magnetic flux (b), the poloidal 
magnetic field (c) and the Poynting flux (d). Each curve describes the structure of the 
core at some distance z. The curves are labeled by the constant C from the left boundary 
condition (11471) : the less C the larger the corresponding z. The dashed and dotted lines show 
the asymptotics oc X and ^ oc X^, Eqs. (11601) and (11581) . correspondingly. 



-1 1 2 3 4 5 

lg(^/7in) 



Fig. 7. — The flux surfaces at k = 2; /? > 1/4; 7max = 1-3 ■ 10^. Tlie dotted line 
shows the boundary of the moderately magnetized core according to Eq. fll22p . h = 
(l/4)lgp-l/4)/«] 




lg(^/7in) + h 



Fig. 8. — Evolution of the Lorentz factor (dashed) and of the ratio of the Poynting to the 
kinetic energy flux (solid) along the flux surface marked by asterisks in Fig. 7. 
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constant C in the left boundary condition fll47p . which could be done by bisection: choosing 
Ci and C2 such that the first solution goes at large to r larger than that of Eq. fl6Tl) 
whereas the second one goes to r smaller than that of Eq. flMl) . one finds the solution for 
C = (Ci + C2)/2 and continues until the necessary solution is found. 

As an example, we presented in Fig. 7 the structure of the jet confined by the outer 
pressure V = PZ~'^; (3 > 1/4. In the Poynting dominated region, the shape of such a jet is 
described by Eq. ( 1551) . One sees that initially the flux surfaces within the jet div erge as Z^^"^ 



together with the boun dary of the jet but eventually a cylindrical core is formed (iTomimatsu 



I994J : iBogovalovl I1995I ) . The transition to the cylindrical flux surfaces occurs when the flow 
becomes moderately magnetized, a ~ 1, and the Lorentz factor saturates, see Fig. 8. The 
transition is in fact very slow therefore the flux surfaces \E' ^ \E' become cylindrical only at 
extremely large distances. The further from the axis, the later the saturation is achieved. 
One sees from Fig.8 that while the main body of the jet remains Poynting dominated, a in 
the core remains larger than unity just approaching unity from above, which agrees with the 
general analysis presented above. 



9. Conclusions 

In this paper, we developed an asymptotic theory of relativistic, magnetized jets. The 
study was motivated by the fact that acceleration and coUimation of relativistic MHD out- 
flows occur in a very extended zone far beyond the light cylinder. This is because the Lorentz 
force is nearly compensated by the electric force when the flow speed approaches to the speed 
of light. Because the dominant terms in the the full set of MHD equations nearly cancel each 
other in the far zone, it is difficult to solve them directly even numerically. In this paper, we 
derived asymptotic equations, which describe relativistic, steady state, axisymmetric MHD 
flows in the far zone. These equations could be easily solved numerically because they do not 
contain either intrinsic small scales like Qr or terms that nearly cancel each other. Moreover, 
in many cases one can solve them analytically or semi-analytically and find simple scalings, 
which provide qualitative understanding of the basic properties of relativistic MHD fiows. 

We applied these equations to externally confined, collimated fiows. Qualitative anal- 
ysis shows that there are two regimes of collimation, which we called equilibrium and non- 
equilibrium, correspondingly. In the first regime, the fiow structure at any distance from the 
source is the same as the structure of an appropriate cylindrical fiow. We call this regime 
equilibrium because in this case, the residual of the magnetic hoop stress and the electri- 
cal force is balanced, as in true cylindrical configurations, by the pressure of the poloidal 
magnetic field. In the non-equilibrium regime, the pressure of the poloidal field is negligibly 
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small so that the flow behaves as if it possesses purely azimuthal field. Such a flow could be 
conceived as composed from coaxial magnetic loops. 

The ourflow is in the equilibrium within the parabola r'^Q < z whereas the non- 
equilibrium regime occurs only outside this parabola, i.e. if the jet is not too narrow. Close 
enough to the axis, the flow is always in the cylindrical equilibrium. An interesting feature 
is that even though the pressure of the poloidal field does not hinder the collimation in the 
non-equilibrium regime, collimation is in fact slower in this regime than in the equilibrium 
one. The reason is that one can neglect the poloidal field only if the flow expands rapidly 
enough. In the two collimation regimes, the flow is accelerated in different ways. In the 
equilibrium regime, the flow Lorentz factor goes as 7 ~ f2r whereas in the non-equilibrium 
regime, the scaling is 7 ~ a/T^/t ~ 2;/r. 

The shape of the flux surfaces in the Poynting dominated, externally conflned jet could 
be found by solving a simple ordinary differential equation. We studied in detail the structure 
of jets with a constant angular velocity conflned by the external pressure with the power 
law proflle, p oc z'^^ . At k < 2, the jet acquires a parabolical shape r oc 2;^, where < 1 
depends on the pressure proflle. The jet is collimated and accelerated until the flow ceases 
to be Poynting dominated. The larger the initial a, the larger the flnal Lorentz factor of 
the flow. The opening angle, B, decreases continuously so that the flow remains causally 
connected, O7 < 1. At k > 2, the flow becomes asymptotically radial. If k only slightly 
exceeds 2, the flow still could be collimated before the flow lines become straight. The flnal 
collimation angle depends only on the pressure proflle. When the flow becomes radial, the 
Lorentz factor could continue to grow so that the flow could become causally disconnected, 
G7 > 1. However, the acceleration is practically saturated when flow re aches the terminal 



Lorentz factor 7^ ^ 7maxO '^1'^ . This generalizes the well known result ( Tomimatsu 1994 



Beskin et al.lll998l ) that the non-coUimated flow is accelerated practically only to 7 ~ 7, 



1/3 

max. 



The Poynting flux generally goes to zero at the axis of the flow therefore a a ~ 1 core 
is always presented in the Poynting dominated jets. We have shown that as the flow is 
accelerated, this core expands and the flow lines within the core approach cylinders. At 
K < 2, the core expands until the a ~ 1 region eventually occupies the whole jet whereas at 
K > 2, the main body of the flow remains Poynting dominated up to logarithmically large 
distances. 
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